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Abstract 



There exists a generic relationship between the thermodynamics of a contin- 
uous polymer with a generic self-interaction and the two-point function of an 
interacting field-theory. In addition, the (2iV)-pomt function of the resulting 
field theory is similarly related to a system of N interacting polymers. 

In the present paper, this relation is explored for the special case of a 
polyelectrolyte, characterized by a screened Coulomb pair potential. The 
corresponding field theory can be recast in a particularly simple form, corre- 
sponding to a quantum-mechanical particle, self-interacting via the emission 
and absorption of a massive scalar field. 

This is particularly useful in a perturbative treatment: perturbative ex- 
pansions for the polyelectrolyte can be simply derived from the loop expansion 
for the related field-theoretical two-point function, for which established com- 
putational methods exist. 
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Introduction Relations between Euclidean field theories and random walks or polymers 
are well-known |]-f|]. For a non- interacting continuous polymer, equivalent to a random 
walk, the relation to the propagator of a free field theory is trivial. For the case of a 
continuous polymer with a delta-function pair interaction, corresponding to a self-avoiding 
random walk, de Gennes has shown the relation to a 4 field theory (Landau-Ginsburg 
model) [0-0. A similar relation has been exploited for an unscreened polyelectrolyte J7j. 
The focus of the present paper is to exploit the generic result for the special case of a screened 
polyelectrolyte, leading to a specific, very simple field-theory. 

The generic interacting polymer A quite generic self-interacting discrete polymer with 
N monomers in .D-dimensional space is described by the Hamiltonian 

1 ^-1 N 

H (xi, . . . , xjv) = - ( x *+i ~ x *) 2 + J2 Vl ( x ^) + 
1 i=i i=i 

N N N N N 

J2 J2 V* (Xi, Xj) + J2 J2 ^3 (Xi, Xj, Xfc) + . . . , 

i=l j=i+l i=l j=i+l k=j+l 

where Xj denotes the position in TZ D of the ith monomer. The first term represents a 
harmonic bonding force between nearest-neighbour monomers, while V\ (x) is a generic one- 
particle potential, V<>(x, y) a ditto pair potential, etc. The interaction is the most gen- 
eral one that depends on positions only and is independent of monomer indices. Ther- 
modynamic averages at temperature T are defined based on the Boltzmann distribution, 
P (xi, . . . , x^ oc exp (-H/T). 

For the Boltzmann exponent S = H/T, we obtain in the continuum approximation, 
defined by replacing the discrete monomer index i by a continuous chain parameter t = Ti, 
such that Xj — > x(Tz) = x(r), the following expression: 

S[x(t)] = 1 Tdr± 2 + ["dr V x (x(r)) + 
2 Jo Jo 

\ dr /Vv 2 (x(t),x(t')) + 
2 Jo Jo 

\ r dr r d^ r ^ v 3 ^M^Mr")) +..., 

b Jo Jo Jo 

where square brackets [] are used to indicate that S is a functional. The length of the 
polymer is represented by /i = NT, while = Vk/T k+1 , for k = 1,2, . . .. The dot in x 
stands for the derivative with respect to r. The continuum approximation is valid when the 
nearest-neighbour distances are small in comparison to the range of the interactions. The 
independence on monomer index is turned into a translation invariance with respect to the 
chain position r: V n depends on (r, r', . . .) only via (x(r), x(r'), . . .). 

The continuum formulation will be the starting point for the analysis carried out in the 
rest of this letter, and 5[x(r)] will be referred to as the action of the system. To begin with, 
we will assume fixed B.C., with the endpoints of the chain locked at x(0) = a, x(/i) = b. 
The main object of interest is the partition function, which is given by the functional integral 



Z(p,a,b) = f (M) b rf[x(r)]exp(-S[x(r)]). 

Jx(0)=a 



The proper partition function for free boundary conditions is obtained simply by integrating 
Z(fi, a, b) over the endpoints a, b. 

The relation to a field theory The generic relation to a field-theory is has been described 
in the literature (see e.g. || and references therein). For the sake of completeness, we will 
here provide a derivation. 

The translation-invariance of the interaction ensures that the interaction part S[ n i of the 
action be a functional only of the particle density, p(y) = Jq dr8 (y — x(r)) , in terms of 
which it is given by 



dxp(x)Vx(x) + -J dxp(x) J dyp(y)^ 2 (x,y) 
- J rfxp(x) J rfyp(y) J dzp(z)y 3 (x, y, z) + . . . 



6 

as a functional Taylor expansion. Next, we assume that the corresponding Boltzmann factor, 
exp(— Si n t), can be functionally Laplace transformed with respect to the (non- negative) 
particle density, 

F[£(x)] = J d[p(x)]exp {-J dx£(x)p(x) - S int [p] 

This is well-defined for any combination of arbitrary one-particle potentials and non-negative 
many-particle potentials. The functional Laplace transform can be inverted to yield, for- 
mally, 

exp(-Si nt [p]) oc j d[f(x)]F[f(x)]exp (J dx£(x)p(x] 

where each ^ (x) is to be integrated along a complex contour to the right of all singularities. 
We then obtain for the partition function 



Z(p,a,b)cx J d[£(x)]F[£(x)] x 

/ d x r exp -- / dr x 2 + / dr£(x(T , 

A(0)=a V 2 Jo Jo J 



and — £(x) appears as an effective one-particle potential for the polymer. Now the x in- 
tegral is recognized as the Euclidean quantum-mechanical path integral for the amplitude 
< a| exp(— piT)|b > to go from a to b in imaginary time p, with the quantum-mechanical 
Hamiltonian 



H = -^V 2 - f (x) + const. 



We thus have 



Z(//,a,b) oc J d[£(x)]F[£(x)] (a |exp 
Now, we can make a Laplace transform in /x, 

/•OO 

W(s, a, b)= / d//exp(— s/i)Z(/i, a, b), 

JO 

yielding 

W(s,a,b) oc / d[£(x)]F[£(x)] (a (if + s)^ 



The matrix-element can be written as a functional integral over a field 0(x) having n = 
components, 



and we obtain 



2 J d[0(x)]0i(a)0i(b) x exp ^- ^ dx0 (# + s) 

W(a,a,b) oc | d[£(x)]F[£(x)] / d[£(x)]&(a)&(b) x 
exp (- 1 dx |i (V0) 2 + s$ 2 - £(i 



Now we can do the £ integral to recover exp (— <Si n t[p( x )])i but with a modified particle 
density, p(x) — > 0(x) 2 . Thus, we finally arrive at the following result: 

W(s,a,b) oc J dtfWM^M*) exp (-5[0(x)]) . 

The RHS defines the two-point function G(s, a, b) of an interacting field theory, with an 
action S given by 

5[0(x)] = |dx{i (V0) 2 + S 2 (x) + Vi(x)^(x)} + 
Ifdxf dy V 2 (x,y)0 2 (x)0 2 (y) 
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+ 
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dxj dyj dz \/ 3 (x,y,z)0 2 (x)^(y)0 2 (z) + . . . 

This shows that the partition function for a generic interacting polymer, when Laplace 
transformed with respect to the chain-length /i, yields the propagator of a zero-component 
field theory, with a bare mass determined by the corresponding conjugate variable s, and 
where the interaction appears in terms of a particle density given by the square of the field. 
Upon collecting the neglected constants of proportionality accumulated in the tranforma- 
tions, the exact relation is given by 

Z(/a, a, b) 



2G(s, a,b) 



d/xexp(— s/i)- 



JdaZ (/j,, a,0) 

where Z denotes the partition function of a non-interacting polymer; the correctness can 
be verified by considering the case of no interactions, since the constants of proportionality 
involved do not depend on the interaction. 



Several polymers Having established the relation between the partition function of an 
interacting polymer and the two-point function of a related field theory, one might wonder 
what is the relevance of the higher order moments in <fi. It turns out (the proof is analogous 
to the one-polymer case) that the four-point function 

(0i(a)^(b)^(c)&(d)) 

is proportional to the Laplace transform of the partition function for two polymers, one with 
its endpoints in a and b, the other in c and d, interacting both internally and with each 
other, and allowed to exchange monomers with each other: The parameter /x, conjugate to 
s, merely fixes their total length. In addition, the moment 

corresponds to the sum of the three possible endpoint rearrangements of the two polymers. 
Similarly, the six-point functions are related to a system of three polymers, etc. 

In contrast, a system of N polymers of fixed lengths \ii corresponds to, upon Laplace 
transforming each length individually, a field theory with distinct squared masses 2s, for iV 
subsets (of zero components each) of the field components. 

For dominantly repulsive interactions, the consideration of several polymers will of course 
make sense only in a finite volume (or finite polymer concentration). 

The screened poly electrolyte In the physically interesting case of a screened polyelec- 
trolyte in three dimensions having only a pair interaction given by a screened Coulomb 
potential, 

V^x, y) = g exp(— Kr)/4irr, with r = |x — y|, 
we have the resulting action 

5[0(x)] = |dx{i(V0) 2 + 



s</> 2 (x) > 4 



9 fj fj exp -Kx-y ) - 2 

7T / rfx / d y — n 1 — & ( x )^ (y)- 



Noting that 



2 J J 47r|x-y| 



exp( ~; /c|x " y|) = (-v 2 +« 2 )" 1 (x, y ) 

47r|x — y| v ' 



we can make a further simplification by introducing an additional scalar field, the exchange 
of which mediates the interaction. This gives 



W(s,a,b) cx J d[0(x)] J d[V>(x)]^i(a)0i(b) x 
exp (-S[0(x),^(x)]^ 



with a local combined action S given by 



S[0(x),V(x)] = 

J rfx 1 1 (V0) 2 + S 2 + i (V^) 2 + i«V 2 + ^<?V 

— * 

The variable s, conjugate to the chain length /i, appears only in the (bare) mass, while 
the inverse range k of the interaction potential appears as the ip mass; for the case of an 
unscreened potential, the -0-field will be massless. 

Now we can modify the boundary conditions; noting that the polymer is translationally 
and rotationally invariant, the partition function Z(/i, a, b) can depend only on R = |a — b|. 
Integrating also over R would yield the relevant partition function for the case of free 
boundary conditions. However, it adds no complication to include a factor exp(iP ■ R) in 
the integral, yielding 

/roo 
dR J dfi exp(-s/i + iP ■ R)Z(fi, 0, R) 

roc 

= / dfj,exp(—s/j,)Z(n,P) 
Jo 

oc J dRexp(iP ■ R)G(0, R) = G(P). 

— * 

This yields the two-point function in momentum space, at momentum P, for the field. 
It contains all relevant information on the distribution of the end-to-end distance R of the 
polymer. Thus, eg., the average end-to-end distance squared, (R 2 ), can be obtained from 

Z(v,P) P 2 (R 2 ) , 



Z((i,0) 



= 1 ^ — '- + C(P 4 



Relation in terms of perturbation theory A perturbative (in the coupling g) expansion 
for the polyelectrolyte is obtained by applying an inverse Laplace transform to the loop 
expansion for the resulting local field theory, the latter amounting to the well-known com- 
putation of Feynman graphs. In three dimensions, no ultraviolet divergences will appear, as 
can be seen from a simple dimensional analysis. 

The Feynman rules for computing G(P) are simple: On a single line add lines of 
emitted and re-absorbed ip lines joining the first line in vertices. A <fi propagator carries 
a factor l/(q 2 + 2s) and a -0 propagator a factor l/(k 2 + k 2 ). For every vertex add a 
factor —2i^fg. The external momentum is given by P, and momentum is conserved at 
every vertex. For every loop, add a momentum integral / d 3 p/(27r) 3 . It is easy to see from 
dimensional counting, that in three dimensions the result will be ultra-violet finite for all 
graphs. Note that no 0- loops are allowed (due to n — 0); this can be interpreted as the 
0-field corresponding to a quantum-mechanical particle that is not second-quantized; thus, 
there can be no creation or annihilation of the corresponding particles. 

Thus, to lowest order G(P) is given by the bare propagator 



which is indeed half of the Laplace transform of exp(— /zP 2 /2), reflecting the distribution of 
R for a free polymer being proportional to exp(— R 2 /2/i). To first order in g, corresponding 
to one loop, there is only one graph contributing, displayed in Fig. [I], yielding 



G 



-Ag f d 3 q 1 



(P 2 + 2s) 2 J (2tt) 3 (q 2 + k 2 )((P - q) 2 + 2s) 

-Ag 1 ( P \ 

■ arctan 



(P 2 + 2s) 2 4ttP \^Ts + k)' 

To second order there are three distinct graphs, displayed in Fig. |2|. 

To compute the partition function Z, and related entities like {S- m \), it is enough to 
consider P = 0, while for (P 2 ), the leading small P 2 corrections must be taken into account. 
Work on perturbative calculations for the polyelectrolyte relying on the above analysis is in 
progress, and will be presented elsewhere [T(| . 



In a similar way, any pair-potential that only depends on the relative position, V^x, x') = 
U (x— x'), can be perturbatively treated; provided the if) propagator is replaced by the Fourier 
transform of the potential, Z7(q) = / <ixexp(iq • x)P(x). 
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FIGURES 




FIG. 2. The three graphs contributing to second order in the perturbation expansion for G(P). 



